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We analyze the effect of typical, unknown perturbations on the 2D toric code when acting as a 
quantum memory, incorporating the effects of error correction on read-out. By transforming the 
system into a ID transverse Ising model undergoing an instantaneous quench, and making extensive 
use of Lieb-Robinson bounds, we prove that for a large class of perturbations, the survival time 
of stored information grows at least logarithmically with the system size. A uniform magnetic 
field saturates this scaling behavior. We show that randomizing the stabilizer strengths gives a 
polynomial survival time with a degree that depends on the strength of the perturbation. 
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The theories of quantum error correction and fault tol- 
erance prove that quantum information can be stored 
despite the deleterious effects of external noise and ex- 
perimental imperfection ^. However, the massive over- 
heads induce significant practical obstacles [2|. In con- 
trast, classical memories are stable without active error 
correction. Can quantum information also be passively 
protected? To this end, quantum memories [31 [3], in 
which a qubit is encoded in the degenerate ground space 
of a Hamiltonian, have generated significant interest re- 
cently. One should design the Hamiltonian structure to 
prevent catastrophic accumulation of errors such that, 
after storage, a single round of error correction correctly 
returns the initial state. Crucially, the aim is to achieve 
storage times that scale with the system size. 

The toric code of 2N^ qubits in 2D [S g] is the proto- 
typical proposal of such a system. While there is strong 
evidence that it is not a good memory at any finite tem- 
perature E] , studies can yield significant insights for 
other types of noise, such as the imperfect implementa- 
tion of an experiment at zero temperature, i.e. the effects 
of unknown static perturbations to the toric code Hamil- 
tonian. The topological phase of the toric code is robust 
against perturbations [3 [3 |H] , meaning the degeneracy 
is only lifted to an exponentially small degree (in TV), 
provided the perturbation strength (5 ^ A, the unper- 
turbed gap. States stored in this ground state space take 
exponentially long to dephase. If the perturbation is un- 
known, the challenge is to encode in this space, which 
may be achieved via an adiabatic path |9l[T0]. At best, 
the final state is subject to a finite density of anyonic ex- 
citations which perturbations can easily propagate into 
logical errors [11] |T2] , although randomizing the weights 
in the unperturbed Hamiltonian induces Anderson local- 
ization and should reduce propagation in almost all cases 
[13] , while leaving the worst-case scaling unchanged [11] . 

Non-adiabatic methods, assisted by error correction, 
focus on accurately preparing the toric code. If this is 
not the ground state of the perturbed Hamiltonian, after 
some time the evolution may mask the stored informa- 
tion. In [HI [12] pathological local perturbations {S <^ A) 
showed the worst case survival time is no better than 



0{5^^ log(A^)). This paper considers a more typical ex- 
perimental affliction, that of a uniform magnetic field. 
This is achieved by transforming the model into a parallel 
set of ID transverse Ising chains subject to an instanta- 
neous quench. While this model has previously been an- 
alyzed in the thermodynamic limit |14| , we study the be- 
havior for large, but finite, system sizes TV, showing that 
the stored data is stable for times ~ A^S~^ log N. We 
also analyze a system with randomized strengths, prov- 
ing a polynomial survival time. 

The Toric Code is defined for an N x N periodic square 
lattice with a qubit placed in the middle of every edge, as 
depicted in Fig. [T] The Hamiltonian is a sum of 4-body 
commuting terms, [Kn,K„i] = 0, 



H=-Y,\^Kn■ 



These terms Kn are typically ZZZZ on a face or XX XX 
around a vertex, where X and Z are the standard Pauli 
matrices. There is a 4-fold degeneracy in the ground state 
space (defined by Kn \ip) = \fp)), allowing the encod- 
ing of two qubits. We consider a rotated version (apply 
Hadamards along every second row) in which every stabi- 
lizer term is the same, XZZX [53]. In addition to mak- 
ing the model translationally invariant (if A„ = A > 0) , 
the logical Z (X) rotations for the encoded qubits are 
two inequivalent columns (rows) of Z (X) operators. The 
Hamiltonian is subsequently subject to a perturbation 
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Our aim is to determine times and field strengths such 
that, with high probability, error correction on the state 
^-it{H+v) 1^^^ produces a logically X rotated state, IV'i) 
{h — 1). This is the time at which, for that model of 
error correction, the stored data is unreliable. 

To locate errors, we measure the stabilizers, and per- 
form minimum weight perfect matching (MWPM) to de- 
termine how to eliminate them. We assume that this 
process is perfect, otherwise the effectiveness, and hence 
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FIG. 1: The 4-body stabilizers and logical X-rotations of the 
toric code. Anyons are propagated horizontally by X rota- 
tions and vertically by Z rotations. 

storage time, is reduced. Since MWPM is difficult to an- 
alyze, we make use of a special property of the chosen 
perturbation. Specifically, a given X„ commutes with all 
the stabilizers, except the two positioned to the left and 
right of the site n. Consequently, all the excitations of 
the original code due to the perturbation manifest them- 
selves along horizontal lines. Every row must have an 
even number of errors. The chance of this happening by 
accident is vanishingly small, implying that the pertur- 
bation is just creating errors along the rows. Hence, we 
perform error correction by using MWPM independently 
along each row. When suitable, this method is strictly 
stronger than the full 2D error correction. If each row 
has a logical error with probability p, a logical error arises 
overall with a probability i(l — (1 — 2p)^). For large N , 
any finite p is destructive. 

The transverse Ising model: The row-like feature that 
the perturbation induces means that the spectrum of the 
Hamiltonian, along each row, is equal to that of the trans- 
verse Ising model with periodic boundary conditions |16| 
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A term Xn commutes with all ZZ terms other than 
Z„_iZ„ and Z„Z„+i. The bit-ffip distribution along a 
particular row of q-^*-(h+v) |^^^ relative to |-0o) is the 
same as that on e~'*^^ |0) . The equivalent error cor- 
rection measures the stabilizers Z„Z„+i, and minimizes 
the number of flips to return the system to a -1-1 eigen- 
state of the stabilizers. Error correction fails if there is 
a finite probability that more than half the qubits have 
ffipped, i.e. if m, the value of the magnetization 

n=l 

after projection onto the stabilizer space, is negative. We 
wish to argue that certain values of 



correspond to finite p. Lieb-Robinson bounds TT] show 
that at time t correlations between terms separated by 
a distance \i — j\ >2vt are exponentially small (hence 
negligible), i.e. for any operators A and B separated by 
distance cIab, there exist positive constants C, 77, 

\\[e-'"'Ae'"*,B]\\ < C\\A\\\\B\\e'^'^'''~'^^^\ 

where z; ~ J = max Sn is the speed of sound for the sys- 
tem. In this regime, (ZiZj) « {Zi){Zj) [19,, so for times 
t <C N/6, almost all the variables Zi are independent. 
Hence, Hoeffding's inequality applies to a good approxi- 
mation j24] and, consequently, the probability that m < 
is finite iff (M) -- 0(l/%/]V). This is the signature for 
failure of the memory. 

An upper bound: Since X^ = X'^^ is a conserved 
quantity of Hj, the initial state decomposes in terms of 

|Gi/Z±> = (|0)^'^±|l)^'^)/^. 

Given that {ZmXi/} ~ for all n, 

(M) = i ^ Re {GHZ^\e'"''Zne'""'' \GHZ+) . 

n 

Bounding this by the absolute value and using \GHZ_) = 
Zn \GHZ+) reveals that 



(M) < 1 



N / . 



\{GHZ^ 



\e'^""'^"'e-'"''\GHZ. 



which is restricted to the +1 eigenspace of Xl. In terms 
of the Majorana fermions 



C2n-1 
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Hj is bilinear, with terms 2f^nmCnCm, where h is a 2N x 
27V matrix (|2iV + l) = |1)) 



E 



(5„(|2n-l) (2n| - |2n) (2n - 1|) 



N 



- J2 ^n (|2n) (2n + 1| - |2n + 1) (2n|) , 

n=l 

given that evolution is only in the +1 eigenspace of X^. 
We use /i*^°^ to denote the instance of h with Sn ~ 
and A„ = 1. The GHZ state is a projection onto the 
simultaneous eigenspace of fermion pairs. 
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which is a fermionic Gaussian state with covariance 
matrix /i^^) [jg]^ \GHZ+){GHZ+\ = ^(/iW). A state 
Pi = uj{M) evolves under a Hamiltonian H = 



2 Z-^nj, 
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Furthermore, for two states pi — w(Afi) and p2 = uj{M2) 



N=50, 5=0.1 



Tr(pip2) = ^sJdetiMi) det(Af2 + A/f i). 
Applying these rules leads to the conclusion that 

(M) < i^|det(iS„-il)|S 



(1) 



B„ 



„2/itg-2/i„t^(0)g2/i„tg-2ht^(0) 



i/f 



where /i„ = /i((5„ i— )■ — 5„) results from ZnHjZn 
2SnXn being bilinear in fermions. 

Analysis: Eq. (flj) gives an upper bound - if it evaluates 
to 1/\/N, the magnetization cannot be larger, and the 
memory is unreliable. To see that the bound is tight, 
consider the evolution of \GHZ±). Initially, these states 
remain close to the original states, with a low density of 
bit flips described by Q. We need to examine the phase 



= Arg 



f{GHZ+\ Qte-'^/t \GHZ+) 
V(GiJZ_| Qte-'-ff/t \GHZ^] 



If 6* = for all populated Q, then the bound is exact, so 
we want to argue that 9 remains small. Now consider the 
evolution of the initial state |0) . In a given Q-sector 
(assumed to be of low density) , the state becomes 



Q{\GHZ^ 



-'^|Gi7Z_))/x/2. 



(2) 



As 6 increases towards tt/2, it becomes a GHZ state. 
However, Lieb-Robinson bounds show [191 that a ID lo- 
cal Hamiltonian with speed of sound v cannot create an 
A^-qubit GHZ state in a time less than 0{N/v). Thus, 
our bound on (A/) reveals not only when the memory is 
certainly not stable, but also when it is stable. 

Returning to the calculation of det(i(i?„ — 1)), since 
Bn is unitary and -B„ = i?*, the eigenvalue pairs take the 
form e**^' (with the possible exception of e*^ = ±1). So, 

det(i(i?„ -1)) = ^ n(e''" - 1)(^"''" - 1) 

i 

= 2T^]^(l-cos6'i). 

i 

Provided |a;| < 1, |(1— x) < e^^^+^^Z^, which reveals that 
W < ^ E„ exp (=^(Tr(B„) + 2A)). We wiU take the 
trace of Bn using a basis \9^) — (|2m — 1) ± i \2m))/y/2, 
which is a diagonal basis of h'^^\ Lieb-Robinson bounds 
again allow us to describe the maximum distance that 
the \6^) can be propagated by h for any model {A„, (5„}. 
Indeed, provided \m — n\ ^ St, 



-2ih„t 2iht 



i) ~ l^m) 



because the propagating states never 'see' the different 
coupling strength between /i„ and ft,, so this means that 
only 0{6t) of the \9^ may not result in —1 for the trace. 
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FIG. 2: Time evolution of the upper bound to magnetization 
for several values of system size and perturbation strength - 
A„ = 1 and 5n is the same at every site. 

For those states 16*^) whose horizon contains the sites 
2n— 1 and 2n, (0+ | i?„ |0+) simply measures the weight 
of the state g-zi/i-.tgai/it ^q+^ ^^ ^^i^ subspace 16*^). Mov- 
ing to the interaction picture with respect to /i, the only 
remaining evolution is ft,„ — h. Consider the eigenspace of 
h. If Sn = 0, we just have the energies {±A„}, and the 
eigenstates are categorized into two bands by whether 
they are of positive or negative energy (corresponding to 
having support on one of the two subspaces 19 ^, )). By 
adding in the {Sn}, there is additional coupling, both 
inter- and intra-band. Within the band, we completely 
solve for the new eigenvectors, treating inter-band cou- 
pling as a perturbation (A,„i„ 3> 6). The positive band 
has energies in the range [Amin — S, A^ax + S] . The first 
order correction to the eigenvectors in a perturbative ex- 
pansion is 0((5/Amin). Whilst hn — h couples between 
all of these eigenvectors, for times t ^ ^n-ln ^^^ rotat- 
ing wave approximation averages away the coupling be- 
tween the two bands. Hence, if a state has support on 
one particular band, it remains in that band. We have 
just argued however, that these bands correspond to be- 
ing in the two different spaces \9^) and \9^), except for 
0((5/A,„in) corrections to the amplitude. Hence no more 
than 0(((5/Anun)^) weight is lost to the other subspace 
during the evolution. Thus, 

Tr(B„) < -(2A^ - /3St) - (/3 - a(J/A^i„)')5i 

where a, (3 represent undetermined constants from our 
perturbation argument and Lieb-Robinson bound. If 
a(5'^t/A^ijj ^ 1, then we conclude that no model will 
ever have a value of {M) that scales worse than 



-a5^t/(16A^,i„) 



1 - a,53^/(16AiiJ 



For sufficiently weak perturbations, the magnetiza- 
tion drops linearly with time. For stronger perturba- 
tions and sufficiently long times (which always exist if 
a((5/Ainin)^A^ > 16), the magnetization decays exponen- 
tially towards zero, and the information is unreliable once 
it drops below l/\/N. Hence, the system presents a sur- 
vival time of r2(^~^A^j,jlog A^). Every model {A„,5„} 



characterized by (5/Ainin certainly survives this long. Not 
all models propagate information in a way that saturates 
Lieb-Robinson bounds. However, the uniform field case 
A„ = l,6n = S has a finite group velocity and is hence 
expected to exhibit this scaling. Numerical confirmation 
can be seen in Fig. [2] Although this is a relatively weak 
effect at small system sizes (cJ/Amin = 0.2 has a critical 
system size of about N = 600, at which point survival 
is for a time of approximately 1500A~?^|^), the scaling is 
so weak that it could easily dominate other timescales as 
one approaches mesoscopic systems. 

Random Systems: In random systems, the propagation 
of information is suppressed compared to that predicted 
by Lieb-Robinson bounds, enabling much longer storage 
times. To this end, we aim to induce Anderson local- 
ization by randomizing cither the A„ (by designing the 
system) or the (5„ (as might arise naturally), concentrat- 
ing on the former and maintaining the assumptions that 
A„ > and Amin ^ S. Fig. [3] already confirms such ex- 
pectations numerically. Localization embodies two con- 
cepts: spectral localization (SL) requires that every eigen- 
vector I Afc) of h has exponentially decaying tails, i.e. there 
exist positive constants Ck , rjk independent of N such 
that I (n|Afc) | < Cfce~'"=l"~'^l [25], whereas dynamical lo- 
calization (DL) requires that | (m| e~ \n) \ is exponen- 
tially bounded in \m — n\. While DL, which implies SL, 
holds for almost all instances of the XX model with ran- 
dom magnetic field [2U], to our knowledge, neither has 
been proven for the transverse Ising model. However, 
the interaction within each band of our model is iden- 
tical to the XX model. Hence, one can readily show 
that the perturbative inter-band coupling leaves SL in- 
tact. [5T] showed that if DL can be proven for the matrix 
/i, a Lieb-Robinson bound with finite horizon follows for 
H. The same argument applied to SL shows that there 
exist positive coefficients c, C" and rj' such that 



A=l, S random 



[e~'"'Ae'"\B]\\ < min(c|i|, C")iV||A||||B|| 



-v'dAl 



Hence, the system exhibits a logarithmic light cone akin 
to that of [22] . Replacing all three applications of Lieb- 
Robinson bounds in our previous analysis, we conclude 
that sufficiently large systems are stable for a time at 



least t ' 



i^V^^^iiin^ ) for almost all instances. 



Conclusions: When subject to an unknown perturba- 
tion along a particular field direction {X or Z), then, 
in every instance, information stored in the degener- 
ate ground space of the toric code is stable for a time 
^^(^"■^A^jjjlogiV), comparable to the worst-case scaling 
for any perturbation constructed in [11) . The analysis 
has been confirmed numerically for system sizes of up to 
2 X 10^ qubits. We have also examined randomized stabi- 
lizer strengths [T^ , improving upon previous analyses by 
accurately encapsulating the effects of particle creation 
in the perturbation, and of error correction. This has 
enabled us to show that the survival time almost always 
scales at least polynomially in N for finite system sizes. 




FIG. 3: Comparison of 3 cases (iV = 100): A = 1,5 = 0.1; 
5 = 0.1, A„ e [0.5,1]; A = 1, 5„ e [-0.1,0.1]. Random cases 
are sampled uniformly and averaged over 100 instances. 

which is a consideration entirely absent from |13j , and yet 
is of central importance. In the future, proving dynam- 
ical localization of the ID transverse Ising model could 
show that storage is stable for arbitrary times. 

While we have argued that the use of ID error correc- 
tion is well justified, how would this compare to the full 
2D error correction? If the error density remains low, 
typical occurrences will be single localized errors. Where 
the ID error correction has a threshold probability of 
50% for these errors, we have numerically estimated a 
threshold for the 2D error correction to be 22% (as com- 
pared to 11% H] for arbitrary local errors), removing all 
iV dependence from our proof, and reducing our lower 
bounds on the survival time to a constant. 
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